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Abstract 

Objective: Ou\ purpose was to validate an early enhancement tinne point for accurately measuring the myocardial contrast 
partition coefficient (lambda) using dynamic-equilibrium magnetic resonance imaging. 

Materials and Methods: Ihe pre- and post-contrast longitudinal relaxation rates (reciprocal of Tl) of the interventricular 
septum (Rim) and blood pool (Rib) were obtained from fifteen healthy volunteers and three diabetic patients with 
hypertension using two optimized Tl mapping sequences (modified Look-Locker inversion recovery) on a 3-Tesla magnetic 
resonance scanner. Reference lambda values were calculated as the slope of the regression line of Rl^ versus Rlt at 
dynamic equilibrium (multi-point regression method). The simplified pre-/post-enhancement two-acquisition method (two- 
point method) was used to calculate lambda by relating the change in Rl^ and Rlt, using different protocols according to 
the acquisition stage of the post-enhancement data point. The agreement with the referential method was tested by 
calculating Pearson's correlation coefficient and the intra-class correlation coefficient. 

Results:lhe lambda values measured by the two-point method increased (from 0.479±0.041 to 0.534±0.043) over time 
from 6 to 45 minutes after contrast and exhibited good correlation with the reference at each time point (r>0.875, p<0.05). 
The intra-class correlation coefficient on absolute agreement with the reference lambda was 0.946, 0.929 and 0.922 at the 
6'^ 7"" and 8* minutes and dropped from 0.878 to 0.403 from the 9* minute on. 

Conclusions: The time-efficient two-point method at 6-8 minutes after the Gd-DTPA bolus injection exhibited good 
agreement with the multi-point regression method and can be applied for accurate lambda measurement in normal 
myocardium. 
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Introduction 

Many previous studies have shown that the myocardial 
extracellular volume fraction (ECV), also known as the distribution 
volume fraction of contrast (V a) or fibrosis index, is able to provide 
information valuable for characterizing myocardial fibrosis and 
correlates well with histology findings [1-8]. However, under- 
standing of the use of this index in both normal and abnormal 
subjects and its association with a variety of heart diseases needs to 
be improved. A non-invasive, reliable and time-efficient method 
for calculating ECV is highly desirable. The best known strategy 
for calculating ECV is based on estimating the myocardial 
partition coefficient (X) of extracellular contrast agents as follows: 



ECV = X ■ (1 - hematocrit) [3-5,9,10], where the hematocrit can 
be easily determined by hematological examination. 

The dynamic-equilibrium cardiovascular magnetic resonance 
(MR) imaging method has been widely used to calculate X and 
ECV [4,7,10-13] based on the longitudinal relaxation rates 
(Rl = 1/Tl) of the myocardium (Rl„, = 1/Tl J and blood (Rib 
= 1/Tlb) measured both before contrast injection and at plasma- 
tissue dynamic equilibrium of the contrast agent (the two-point 
method). This method is derived from a regression method 
requiring multiple measurements of Rl^ and Rli, at dynamic 
equilibrium after the contrast is administered [2,14] and is more 
time-efficient in the clinical setting. However, the accuracy and 
reliability of estimates using the two-point method compared with 
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those using multi-point regression method have not yet been 
studied. In addition, ECVs calculated using the two-point method 
have been found to change over time [7,12,15]. Therefore, the 
most efficient and reliable timing for the post-contrast Rl 
measurement remains uncertain. 

The Modified Look-Locker Inversion Recovery (MOLLI) has 
been proposed for measuring cardiac Tl during a single breath- 
hold [6,16-18] and has been used to calculate X in some reported 
studies [4,11,12,15]. However, previously proposed MOLLI 
protocols have mostly been derived using 1.5-Tesla MR systems 
[19-21]. Although improvements have been made using 3-Tesla 
systems in some studies [13,22], these protocols may still have 
limitations due to longer tissue Tl at 3 Tesla, including the 
insufficient longitudinal magnetization recovery between the 
inversion pulses, as well as the heart rate (HR)-dependent accuracy 
[13,22]. 

In this study, we first optimized the MOLLI schemes for Tl 
mapping at 3 Tesla using a tube phantom study and then vahdated 
a two-point method MR imaging protocol for accurate and time- 
efficient measurement of X from the dynamic-equilibrium at 3 
Tesla. The potential mechanisms affecting time-dependent X and 
the ECV estimation when using the two-point method were 
investigated using the normal multi-point regression method as the 
reference. 

Materials and Methods 

Ethics statement 

The study protocol was approved by Ethical Review Board for 
Clinical Research of Zhongda Hospital, Southeast University, and 
the study was performed in accordance with the Declaration of 
Helsinki. AU participants gave written informed consent. 

Description of MOLLI schemes 

The MOLLI sampling algorithm is expressed as readouts 
between inversion, with the wait beats in parentheses, such as 
MOLLI3(3)3(3)5 [18]. The inversion times (TIs) in one MOLLI 
scheme were determined by the minimum inversion time 
(Tlmmimum)) invcrslou time increment (TIi„,.remeiit) and incUvidual 
cardiac cycle as described by Messroghh [19]. 

Phantom study 

To test the T 1 mapping accuracy of different MOLLI schemes, 
nine Gd-DTPA agarose samples with T 1 values of approximately 
200 to 2000 msec and T2 values of approximately 170 msec 
(which are similar to the values reported for post-contrast blood 
[15]) were studied. Three groups of MOLLI approaches with 
gradually changing sampling algorithms (Table 1) were tested. The 
reference Tl values for each sample were determined using 
standard inversion recovery-spin echo pulse sequences (IR-SE) 
[23]. A "HR" produced by a simulated ECG signal was set to 70 
beats per minute (bpm). The MOLLI and IR-SE sequences were 
repeated alternately to assess the precision of the measurements. 
Two MOLLI schemes were selected from this section and, along 
with four protocols selected from the literature (i.e., 3(3)3(3)5, 
3(3)5, 5(3)3, 4(1)3(1)2 [13,19,21,22]), were tested further using 
simulated HRs of 40 to 100 bpm. 

A pipeline system phantom was designed to simulate blood flow 
from a peripheral region into the imaging plane during MOLLI 
acquisition and to study how this "flow-in effect" impacted the Tl 
calculation. The homemade device included a diluted CUSO4 
solution in a bottle as the control sample, and the same solution 
flowed through a long silicone tube twining around the bottle, 
which had a Tl value that was similar to blood (Figure lA). The 



Table 1. MOLLI schemes tested in the phantom study. 
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Note: r/minimum 15 the initial inversion time as realized in the first heartbeat of the 
first Look-Locker experiment; Tlj^crement is the shift Tl increment relative to the 
TIminimum for acquisition in the first heartbeat of the second Look-Locker block; 
/-/ represents the cardiac cycles betw/een the adjacent inversion pulses for 



longitudinal relaxation. 



doi:10.1371/journa!.pone.0093124.t001 

flow velocity of the pipeline was set for a transit time of 
approximately 16 sec, which was slightly longer than the duration 
of the MOLLI schemes (approximately 8 cm/sec). Eighteen 
sampling positions along the silicone tube were simultaneously 
mapped within one image of the mid-vertical section (Figure IB). 
Each sampling position had a different flow volume coming from 
out of the inversion field during MOLLI acquisition. MOLLI 
4(7)3, 6(5)3 and 8(3)3 were performed in this section, with each 
using the same inversion intervals (1 1 cardiac cycles) but a 
different number heart beats for sampling (7, 9 and 1 1 beats, 
respectively). 

Human study 

The study subjects included fifteen healthy volunteers (6 males 
and 9 females, aged 19 to 50 years; mean: 36.3 years) and three 
diabetic patients diagnosed with hypertension (1 male and 2 
females, aged 56 to 69 years; mean: 62.7 years). The subjects were 
recruited between March and May of 2012. The diabetic patients 
were recruited to test the clinical feasibility of the proposed 
method in patients with a high risk for developing heart failure. 
Healthy volunteers underwent hematological examinations and 
twelve-lead electrocardiograms and had their histories taken prior 
to the study to exclude those with other diseases involving the 
heart. Clinical information regarding the three patients was taken 
from their latest medical records no more than 2 weeks before the 
study commenced. None of the subjects had general contraindi- 
cations to cardiac MR imaging or Gd-DTPA. 

All subjects received 0.15 mmol Gd-DTPA/kg body weight 
(MagTievist, Bayer Shering Pharma AG, Berlin, Germany), which 
was infused manually at a rate of approximately 0.5 ml/sec, 
followed by 5 ml of saline at the same rate. A mid-ventricular short 
axis plane was selected for imaging. Tl mapping was performed 
using the optimized MOLLI schemes in accordance with the 
imaging workflow shown in Figure 2. The commonly used 
MOLLI3(3)3(3)5 sequences were also performed before and after 
contrast injection in 7 subjects and were compared with the 
optimized MOLLI schemes in vivo. The total scan times were 
typically 55 to 60 minutes, including patient preparation time. Five 
healthy volunteers (1 male and 4 females, aged 29 to 48 years; 



PLOS ONE I www.plosone.org 



2 



March 2014 | Volume 9 | Issue 3 | e93124 



Time-Efficient Myocardial A, Measurement with MRI 



pump 




Tl [ms] 



1600 
1200 
800 
400 
0 




control 

-4(7)3 
-6(5)3 
-8(3)3 



1600 
1200 
800 
400 



7 9 11 13 
sampling site 



15 17 control 



B 



2000 msec 




O'msec 



4000 



[a.u. 


3000 


>. 




0) 






2000 






"to 




1= 
D) 


1000 


W 















«^ • control 








o sampling site 8 




o sampling site 16 



1000 2000 



3000 4000 
Tl [ms] 



5000 6000 



Figure 1. The pipeline system for testing the influence of flow on T1 mapping. A. The home-made device included the following 
components: a dilute bottled solution of CUSO4 as a control, the same solution flowing through a long silicone tube twining around the bottle, a 
pump providing power to the flow and a bucket outside the scanner for buffer; B. A mid-vertical Tl map of the device generated by modified Look- 
Locker inversion recovery (MOLLI) schemes of 4(7)3; C. MOLLI4(7)3, 6(5)3 and 8(3)3 all underestimated Tl at upstream sampling positions (sites 1-11) 
by a similar degree; D. The signal intensity-inversion time fitting curve deviated significantly (red arrows) at sampling site 8, whereas the skewing 
became very small at the downstream sampling site 16 (blue arrow). 
doi:1 0.1 371 /journal.pone.00931 24.g001 



mean: 36.8 years) participated in a repeated study 12 to 16 weeks 
later to assess the reproducibility of the initial results. 

MR imaging parameters 

AH MR scans were performed using a clinical 3.0-TesIa MR 
scanner (MAGNETOM Verio, Siemens AG Healthcare Sector, 
Erlangen, Germany). A six-channel body matrix coil and a six- 
channel spine matrix coil were used for signal reception. The 
parameters for the IR-SE for Tl measurement were as follows: 
repeat time/ echo time = 10000/8.4 msec, number of averages = 1, 
in-plane resolution = 1.8x2.2 mm, matrix size = 192 x 78 and slice 
thickness = 8 mm. The TIs were 25, 50, 100, 180, 300, 500, 750, 
1100, 1500, 2000, 2700, 3500, 4500, 5800 and 7500 msec. The 
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Figure 2. T1 mapping workflow in the human study. Measure- 
ments were repeated twice before the contrast was administered to 
determine Tl measurement precision, and the mean value was used to 
calculate X. 

doi:10.1371/journal.pone.0093124.g002 



MOLLI schemes are listed in Table 1. The following MOLLI 
parameters were identical between the phantom and human 
studies: balanced steady-state free precession readout, fhp angle of 
35°, parallel acquisition technique acceleration factor =2, 6/8 
partial Fourier k-space sampling, identical spatial resolution and 
matrix size with IR-SE and repeat time/ echo time = 2.5/ 1 .06 msec. 
A three-lead wireless vector-cardiogram was used to detect the R- 
wave that triggered MRI acquisition. AH Tl maps of the heart were 
gated to mid-diastole and collected under an end expiration breath 
hold with patients in the supine position. 

Image analyses and data processing 

IR-SE Tl values were calculated by fitting signal intensities (SI) 
acquired at different TIs to a 2-parameter mathematical model of 
SI = SIo ■ I 1— 2 • exp (— TI/Tl) |, using a home-developed 
MATLAB program (Version R2012a, Math Works, Natick, USA). 
Phantom Tl values of the two selected MOLLI schemes were 
fitted to the IR-SE Tl values using linear regression functions. 
These linear functions were used to correct the measured in vivo 
MOLLI Tl values. 

In each Tl map, regions of interest (ROIs) of 2-3 square 
centimeters were selected from the interventricular septum and the 
left ventricular blood pool. The myocardium-blood border was 
carefully excluded from the ROIs to avoid interference from the 
partial volume effect. We obtained the mean Tl and Rl values 
within the ROIs. In each subject, the blood-tissue contrast 
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equilil)rium state was judged by two independent examiners and 
determined using the deviation of the Rln^ versus Rib from a 
linear relationship [14]. Of the two time points suggested by the 
two examiners, the later value was selected as the time at which 
individual erjuilibrium was established. In each subject, the linear 
regression of the Rl^, and Rli, values, which were acquired with 
optimized MOLLI pre-contrast and at 6 (according to the latest 
equilibrium state establishment time of the 18 subjects) to 45 
minutes post-contrast, was used to generate a slope as the 
reference X; this approach is referred to as the regression 
method in this study. The two-point method was applied to 
calculate a series of X values using different post-contrast data 
points as follows: X = ARlm/ARlj,, where ARl = Rl post-contrast 
— Rl pre-contrast- The apparent X before equilibrium establishment 
was also calculated using the two-point method. All calculations 
for X were performed using Tl values corrected according to the 
functions derived from the phantom study. 

Statistical analyses 

Statistical analyses were performed using SPSS statistical 
software (v. 18, Chicago, IL). Numerical data were reported as 
mean values ± standard deviation, except when otherwise 
specified. The coefficient of variation (CoV) for multiple repeated 
measurements was used to evaluate precision. For statistical 
comparisons, the independent-sample t test, paired-sample t test, 
Pearsc)n's correlation coefficient and Bland-Altman plots were 
applied. The intra-class correlation coefficient (ICC) on absolute 
agreement was applied to assess measurement reproducibility and 
agreement. A^<0.05 was considered statistically significant. 

Results 

Phantom study 

Tl values measured using the different MOLLI schemes 
exhibited good agreement with the IR-SE results of the four 
short-Tl agarose samples (between 213 and 892 msec). Varying 
degrees of Tl underestimation were found in the other five 
samples, which all had longer Tl \'alues (between 1107 and 
1982 msec). In general, MOLLI schemes with shorter inversion 
intervals more profoundly underestimated Tl, and the degree of 
underestimation (in terms of the absolute Tl value) increased as 
the Tl values increased (see Figure SI). After considering accuracy 
and breath-hold durations, the following two MOLLI schemes 
were selected for the in vivo study: 4(7)3 for the pre-contrast 
measurement (Tl = 1000 to 2000 msec) and 4(3)3 for the post- 
contrast measurement (Tl = 200 to 1000 msec) (Figure 3). These 
two MOLLI schemes tended to be more HR-independent than 
schemes reported in the Uterature for their respective applicable 
scopes (Figure 4). The following linear regression equations 
described the relationship between the two MOLLI sequence 
values and IR-SE results and were used to correct Tl 
values measured in vivo: T1ir_se = 1-07 xT1molli4(7)3 ~62.4 
msec (Tl = 1000-2000 msec) and T1ir_se= 1.02xT1molij4(3)3 
-15.8 msec (Tl = 200-1000 msec). The CoVs for repeated 
measurements with IR-SE, MOLLI4(7)3 and MOLLI4(3)3 
(repeated 5, 8 and 8 times, respectively) were all less than 0.01, 
except for the MOLLI4(7)3 sequence with the smallest reference 
Tl (approximately 213 msec), which had a CoV of 0.0101. 

In the pipeline system phantom study, the flow caused a 
noticeable impact on the Tl measurements, as expected 
(Figure IC, ID). Underestimation was more pronounced at 
upstream sampling positions (sites 1 to 8) but became neghgible at 
downstream sampling positions (sites 9 to 18). The three MOLLI 
schemes used in this section exhibited similar degrees of 



underestimation, which suggests that increasing the number of 
sampled heart beats did not mitigate the effects of flow-in when 
calculating Tl. 

Human study 

Tl measurement with MOLLI. When using the optimized 
MOLLI method and corrected with the finear regression 
equations, the pre-contrast Tl„ and Tlj, values of all 18 
participants were 1269.3±61.5 msec and 1 793. 4± 101.4 msec, 
respectively. Between 1 and 45 minutes after contrast injection, the 
Tin, values rose from 336.4±96.1 msec to 718.0±47.5 msec, and 
the mean Tlj, values rose from 14 1.5 ±23.8 msec to 
591.9±53.7 msec (Figure 5). The in vivo comparisons between 
the optimized MOLLI schemes and MOLLI3(3)3(3)5 of the seven 
subjects is shown in Table 2. In results similar to the phantom 
study, MOLLI4(7)3 showed higher Tl^, and Tl|, values than 
MOLLI3(3)3(3)5 before contrast agent injection, whereas 
MOLLI4(3)3 showed smaller Tl„i and Tlj, values than 
MOLLI3(3)3(3)5 between 16 and 20 min after contrast agent 
injection. The breath-hold durations were reduced by 2.8±0.3 sec 
and 6.6±0.8 sec, for MOLLI4(7)3 and MOLLI4(3)3 respectively, 
compared with those of MOLLI3(3)3(3)5. The pre-contrast scan- 
rescan ICCs (95% confidence inter\^al) of MOLLI4(7)3 on Tl„ 
and Tib were 0.956 (0.886-0.983) and 0.937 (0.840-0.976), 
respectively, and were comparable with those of MOLLI3(3)3(3)5, 
which produced ICCs of 0.989 (0.920-0.998) and 0.990 (0.949- 
0.998), respectively. 

Dynamic equilibrium establishment and the regression 
method. The deviation from linearity of Rim versus Rl^, was 
significant during the first and second minutes after contrast 
administration and became milder from the third minute on. The 
deviations were negligible in 6 of the 18 cases at the fourth minute, 
and the deviations were negligible in 16 cases at the fifth minute. 
No significant deviations were found from the sixth minute on, 
indicating that states of dynamic equilibrium were established in 
all 18 cases (at 3, 2 and 5 min after contrast injection in the three 
patients). Interestingly, pre-contrast data points were consistently 
below the regression line in each subject, although the deviations 
were very small. An example of the regression approach is shown 
in Figure 6A. For the 1 8 participants, the r^ value when calculating 
X using pre-contrast and 6-45 minutes post-contrast data points 
was 0.993±0.004 (range: 0.985-0.998), which indicated good 
finearity of Rl^, versus Rlj, values between 6 and 45 minutes 
following contrast; this time period was identified as the dynamic 
equilibrium phase. 

Comparison of X. The mean regression method X value of 
the 18 subjects was 0.468±0.042 (range 0.410 to 0.523), whereas 
the two-point method X value increased over time following 
contrast, from 0.479±0.041 at the 6* minute to 0.534±0.043 at 
the 45* minute (Figure 7). An example of the two-point method 
using different post-enhancement data points (6 min, 16 min and 
35 min) is shown in Figure 6B. Pearson's correlation coefficient (r) 
between the two-point method X and the regression method X 
increased with time while equifibrium was being estabfished before 
the sixth minute; the coefficient remained greater than 0.875 from 
the sixth minute on. The ICCs on absolute agreement of the two 
methods were 0.946, 0.929 and 0.922 at the 6*, 7* and 8* 
minutes, respectively, and dropped over time from 0.878 to 0.403 
from the 9* minute on (Figure 8). The absolute biases (±1.96 
standard deviation) for the two-point method at minutes 6 through 
8 (the early stage) were 0.011 ±0.017, 0.014±0.016 and 
0.014±0.018, respectively, as shown using Bland-Altman plots, 
and diose values increased from 0.019±0.024 to 0.066±0.040 
from the 9'*' minute on (see Figure S2). 
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Figure 3. Comparison of tlie optimized and original modified Loolc-Loclcer inversion recovery (IVIOLLI) scliemes. A. MOLLI3(3)3(3)5, a 
scheme in common use, with a duration as long as 17 cardiac cycles; 8. MOLU4(7)3 for pre-contrast use, with a shorter duration of 14 cardiac cycles 
and a longer inversion interval of 1 1 cardiac cycles; C. MOLLI4{3)3 for post-contrast use, with a shorter duration of 10 cardiac cycles and a longer 
inversion interval of 7 cardiac cycles; D. MOLLI4(7)3 improved the accuracy of the T1 measurement in the five samples with longer T1; E. The 
MOLU4(3)3 T1 values were not significantly different from the MOLLI3(3)3(3)5 11 values (p = 0.132, 0.1 15, 0.254 and 0.672 in ascending order of T1 
values). Both significantly but slightly overestimated the shortest three sample T1 values and underestimated the sample T1 of 892 msec. Heart 
rate = 70 bpm. 

doi:10.1371/journal.pone.0093124.g003 



Reproducibility. In a repeat study of 5 individuals per- 
formed more than 1 2 weeks later, the mean interval separating the 
two studies was 108 days (range: 91 to 119 days). The ICC of X, on 
absolute agreement between the two separate studies was 0.873 
using the regression method and between 0.688 and 0.943 at 
different stages using the two-point method (from 6-45 min). 
Con.sidering a long time interx^al between measurements may 
involve changes in hematocrit, X value.s were individually adjusted 
as follows: ECV — A, • (1 — hematocrit), where the hematocrit was 
measured 0.5-3 hours before each study. The study-restudy ICCs 
for ECV were then calculated, and the values for both the 
regression method and the two-point method were greater than 
0.887 (see Figure S3). 

Discussion 

In this study, we tested and verified the two-point method for 
calculating the partition coefficient at dynamic equilibrium after a 
quick intravenous bolus injection of a contrast agent. With the use 
of MOLLI, X can be calculated within two single breath-hold 



scans before and after introduction of the contrast agent. One 
major finding of this study is that the results of the two-point 
method at the early stage of enhancement (6-8 min after contrast 
agent injection) exhibited excellent agreement with the results of 
the regression method. This finding suggests that properly 
decreasing post enhancement delay time not only saves scanning 
time but also potentially improves accuracy. 

Equilibrium establishment 

Previous studie.s using the two-point method have not generally 
included measurements between minutes 5 and 10, likely to ensure 
fuU dispersion and equilibration [4,6-8,15,24]. Lee et al. reported 
that dynamic equilibrium between the plasma and myocardium 
was reached at 8.5 min [13]. In that study, however, only three 
time points were established during the first 10 minutes (3.5, 5 and 
8.5 min), and equilibrium may have been reached in less than 
8.5 min. The present study sampled more time points during the 
first 10 min following Gd-DTPA injection, and the results 
indicated that equilibrium can be established by the sixth minute; 
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Figure 4. Comparison of different modified Look-Loclter inversion recovery (MOLL!) sampling algorithms with inversion recovery- 
spin echo (IR-SE). y axis: relative T1 values (± SD) to those determined by IR-SE, as measured with different MOLLI schemes on sample T1 of 
1769.9 msec (A), 1249.8 msec (B), 821.6 msec (C) and 218.4 msec (D). The measurements were repeated seven times. MOLLI4(7)3 exhibited the best 
heart rate (HR)-independent accuracy of the six MOLLI schemes on the two longer Tl samples, which were closed to the native blood and 
myocardium Tl. For shorter Tl, all exhibited good HR-independent accuracy except for 4(1)3(1)2. 
doi:l 0.1 371 /journal.pone.00931 24.g004 



this fmdiiig was true in all 18 cases. In the present study, we included 
3 diabetic patients who had peripheral atherosclerosis and risk 
factors for microcirculation dysfunction in the myocardium, and 
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Figure 5. T1 of the myocardium and blood of the 18 subjects 
before and after Gd-DTPA injection. Rings, Tl of the myocardium; 
dots, Tl of blood; bars, mean ± standard deviation. 
doi:10.1371/Journal.pone.0093124.g005 



dynamic equilibrium was also established within sixth minutes in 
these patients. Jerosch-Herold et al. reported that the apparent 
agreed, within 1 % , with the true EC V at times longer than 3 min 
following contrast bolus injection for myocardial blood flow (MBF) 
above 0.5 ml-min '-g ', using a spatially distributed model [14]. 
Even at an MBF of 0.25 ml-min~'-g~', the apparent Vj conformed 
well to the true ECV after 5 min. Compared with most reports, an 
MBF of 0.25 ml-min~'-g~' is extremely low for non-infarcted 
myocardium [25-29], which suggests that an early enhancement 
two-point method might be useful for a wider range of patients than 
only the normal population. In patients with severe ischemic 
myocardium and expected MBF lower than 0.25 ml-min ' -g ' , an 
additional scan at 10 to 15 min would be helpful to ensure 
equilibrium. 

Time dependence 

Consistent with previous studies [7,12], the results of the present 
study indicate that X (or ECV) continuously increased over time 
when using the two-point method. This phenomenon has been 
observed in other studies, but the underlying reasons have not 
been discussed [15,30]. Because the accuracy of MOLLI depends 
on measuring true Tl values and MOLLI tends to underestimate 
longer Tl times, the pre-contrast point is more likely to deviate 
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Table 2. In vivo comparison between optimized MOLLI and MOLLI3{3)3{3)5. 









Reference 


MOLLI4(7)3 


MOLLi4(3)3 


MOLLI3(3)3(3)5 


Duration (sec) 




N.A. 


13.2±1.6 


9.4±1.2t 


16.0±2.0tt 


Inversion pulse interval (sec) 




N.A. 


10.4±1.3 


6.6±0.8t 


5J±0.7tt 


Pre-contrast T1 (msec) 


Myocardium 


1267.1 ±40.1 


1242.5±37.5* 


N.A. 


1212.1±31.1*t 




Blood pool 


1823.7±90.6 


1762.7±84.7* 


N.A. 


1708.6±86.7*t 


Post-contrast T1 (msec) 


Myocardium 


586.7 ±35.4 


N.A. 


590.7±34.7* 


602.8±30.0*t ^^^^ 




Blood pool 


426.5±25.7 


N.A. 


433J±25.2* 


439.5±23.7*t 



Note: Data are given as the mean ± standard deviation, n = 7. Reference T1 values were measured using optimized MOLLI and corrected by the corrected function 
derived from the phantom study. Post-contrast T1 values were calculated as the mean values of 16, 18 and 20 min for MOLLI4(7)3 and the mean values of 17 and 
19 min for MOLLI3(3)3(3)5. Significant differences In the paired sampling f test (p<0.05) are labeled by * when compared with the reference, f compared with the 
MOLLI4(7)3, and t compared with the MOLLI4(3)3. 
doi:10.1371/journal.pone.0093124.t002 



from the theoretical linearity between Rl^, and Rlj, (the primary 
cause of this deviation may be the underestimation of blood Tl, 
which will be discussed later) and leads to a time-dependent, 
biased estimation (Figure 6). This mechanism may explain the 
discrepancy among the time dependence values of X (or ECV) 
reported in studies using different Tl measurement protocols [4]. 
However, the estimated apparent X also increased as dynamic 
equilibrium was being estabhshed, but this increase cannot 
account for persistent time dependence after such a long delay 
time. Earlier post-enhancement acquisition should increase the 
accuracy of X estimation when using the two-point method, as long 
as equilibrium has been established, as exhibited by the better 
agreement between the results of the two-point and multi-point 
regression methods in our experiments. The reason for this 
increased accuracy may be the greater ARl at the earlier 
enhancement stages lowered the impact of the inaccurate 
precontrast Tl measurement. Alternatively, the injection of a 
higher dose of contrast agent is also helpful for increasing ARl and 
therefore improves X estimation accuracy even when post-contrast 
images are acquired much later than equilibrium is established. In 



another word, delay time and contrast dose can impact the 
accuracy of ECV estimates by the same way. This finding may 
partially explain the higher ECV estimates in the 0.1 mmol/kg 
group compared with the 0.15 and 0.2 mmol/kg groups in 
Miller's study, despite the use of a uniform and extended delay 
time [7]. 

Optimized MOLLI schemes 

Previous studies have discussed the optimization of MOLLI 
parameters and their reproducibility using a 1.5-Tesla MR 
scanner [18,19], but those studies recommended MOLLI schemes 
with an inversion interval of 6 cardiac cycles. This interval may be 
insufficient for longitudinal magnetization recovery at 3 Tesla due 
to the increased Tl at higher field strengths [31]. A significant Tl 
underestimation was reported in a previous study by Messroghli, et 
al., caused by the short inversion interval [18]. Moreover, HR- 
dependent inversion intervals change from 4000 to 7200 msec in 
the HR range of 50-90 bpm when these MOLLI schemes are 
used, which makes HR-based Tl correction necessary [13,18]. In 
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Figure 6. Multi-point regression and two-point methods. A. The multi-point regression method at dynamic equilibrium; B. The same data 
were then used to calculate X using the two-point method with representative post-enhancement time points of 6 min, 16 min and 35 min after 
contrast administration. The slope of the straight line determined by the corresponding two points was used to represent the three two-point 
method X values. The deviation from the regression linearity of the pre-contrast point resulted in an increase in the X measurement over time, which 
was observed with the two-point method. Data were taken from a 56-year-old man with diabetes and hypertension. 
doi:10.1371/journal.pone.0093124.g006 
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doi:1 0.1 371 /journal.pone.00931 24.g007 



contrast, much longer inversion intervals of 1 1 and 7 cardiac 
cycles (for MOLLI4(7)3 and MOLLI4(3)3, respectively) were 
selected in our study to allow for the full recovery of longitudinal 
magnetization, and as a result, the accuracy was improved. The 
inversion intervals are more than 7300 and 4600 msec for 
MOLLI4(7)3 and MOLLI4(3)3, respectively (more than 5.5 times 
corresponding Tl^, and post-enhancement Tl|, values at 3T), with 
HR values less than 90 bpm. At that range, the influence of HR 
becomes negligible; however, underestimation may still be an issue 
when measuring pre-enhancement Tli,. To compensate for 
extended sequence duration, we reduced the number of Look- 
Locker inversion recovery blocks to shorten breath-hold duration. 
Despite the corresponding reduction in sampling points, particu- 
larly in the first heartbeat after inversion, the precision was 
comparable with that of the reference IR-SE. 

Half-inverted inflow effect 

We demonstrated, using a pipeline device, that Tl could be 
underestimated because of imperfect inverted blood flow to the 
imaged slice, similar to the inflow effect encountered when using 
dynamic contrast-enhanced MR imaging [32] . The underestima- 
tion of Tl depended on many factors, including flow velocity, flow 
direction, sequence type and repeat time. Nevertheless, we found 
that the Tl values measured at the downstream sampling positions 
(sites 12-18) were very similar to the values measured in the stable 
solution sample (Figure IC), which suggests that the influence of 
other factors was negligible. Although a non-spatially selective 
inversion pulse was used, the volume excited by the inversion pulse 
(approximately 40-50 cm, depending on hardware) may still have 
been limited. During acquisition, half-inverted blood outside the 
inversion volume might flow into the right atrium and right 
ventricle (upstream positions) through the cava, producing 
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Figure 8. Agreement between the two-point method at 
different stages and the regression method. The agreement 
between the two methods was evaluated with an intra-class correlation 
coefficient (black) and Pearson's correlation coefficient (white), n = 18. 
doi:10.1371/journal.pone.0093124.g008 



inhomogeneous Tlj, in the right atrium and right ventricle 
(Figure 9). Blood in the left atrium and left ventricle (downstream 
positions) came directly from the kmg and right heart within 4 
heart beats, which was adequate for inversion volume coverage. It 
is not difficult to infer from the signal recovery fitting curve of fluid 
(Figure ID) that blood flow with longer Tl is more likely to be 
underestimated than that with a shorter Tl. The significance of 
these results lies in the following hypotheses: (a) the blood Tl value 
in the heart is underestimated, (b) that underestimation may 
become insignificant when MOLLI is performed after contrast 
injection and (c) the ROI used to measure blood Tl should be set 
in the left ventricle, where little fresh blood may flow in during 
MOLLI acquisition, whereas the right ventricle is more susceptible 
to imperfect inverted blood flow. 

Reference method 

There is still no gold standard MR protocol for calculating X. A 
bolus of the extracellular contrast agent gadolinium (Gd-DTPA) 
followed by continuous infusion has been recommended in some 
previous studies as a method for achieving blood-myocardial 
contrast equilibrium, with a constant concentration maintained 
during post-contrast scanning [3,8,33]. This enhancement tech- 
nique, however, still might bias the X due to underestimation of the 
pre-contrast Tli,. Equilibrium in our study was defined as a 
dynamic process, with the equilibrium concentration of the 
contrast dropping slowly due to renal excretion. This equilibrium 
allowed X to be calculated with a wider range of Rl^ and Rli, 
values using the regression method, which limited possible 
systemic errors caused by pre-contrast Tl measurement. 



2000 msec 




0 msec 



4-chamber view 



short-axis view 



Figure 9. A test pre-enhancement T1 mapping with MOLLI4(7)3 
using a 4-chamber view (left) and a short-axis view (right). The 

black arrow indicates the heterogeneity of blood Tl in the right atrium 
(RA). The white dashed lines indicate the relative location of each map. 
LV, left ventricle; LA, left atrium; RV, right ventricle. 
doi:10.1371/journal.pone.0093124.g009 
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Study limitations 

This study has some limitations. First, there is still no gold 
standard method for measuring the Tl of blood. Our pipeline 
system demonstrated the effects of inflow on the T 1 measurement, 
and we discussed how to reduce this influence. Second, no patient 
with definite scarring or severe ischemia of the myocardium was 
included, which limits how far our conclusions can be extended. The 
applicability of the two-point method at early enhancement stage 
within different myocardial injuries needs further confirmation. 

Conclusions 

In conclusion, optimized MOLLI4(7)3 and 4(3)3 parameters are 
reliable for clinical use. Measurements taken using a time-efficient, 
two-point MRI technique at 6-8 minutes after a Gd-DTPA bolus 
injection exhibited good agreement with the multi-point regression 
method and can be applied for accurate X and ECV measurements 
in normal myocardium. Tl -dependent MOLLI accuracy is a key 
factor for the time dependence during X and ECV measurements 
using the two-point method at dynamic equilibrium. 

Supporting information 

Figure SI Comparison of the three groups of MOLLI 
schemes with the IR-SE reference. MOLLI exhibited good 
agreement with the IR-SE reference for shorter Tl values and 
underestimated longer Tl values by difF(;r(;nt degrees, which was 
associated with the inversion intervals. Heart rate = 70 bpm. 
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